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We study the properties of tiie cliromoiiydrodynamical instabilities induced by a relativistic jet 
that crosses the quark-gluon plasma. Assuming that the jet of particles and the plasma can be 
described using a hydrodynamical approach, we derive and discuss the dispersion laws for the 
unstable collective modes. In our analysis the chromohydrodynamical equations for the collective 
modes are tackled in the linear response approximation. Such an approximation, valid for short 
' time scales, allows to study in a straightforward way the dependence of the dispersion laws of the 

I collective modes on the velocity of the jet, on the magnitude of the momentum of the collective 

, mode and on the angle between these two quantities. In the conformal limit we find that unstable 

' modes arise for velocity of the jet larger than the speed of the sound of the plasma and only modes 

' with momenta smaller than a certain values are unstable. Moreover, for ultrarelativistic velocities 

of the jet the longitudinal mode becomes stable and the most unstable modes correspond to relative 
angles between the velocity of the jet and momentum of the collective mode larger than ~ 7r/8. 
Our results suggest an alternative mechanism for the description of the jet quenching phenomenon, 
00 \ where the jet crossing the plasma loses energy exciting colored unstable modes. 
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I. INTRODUCTION 



(— I I One of the remarkable findings at the Relativistic Heavy Ion Collider, RHIC, is that hydrodynamical models are 
able to accurately describe the spectra of soft hadrons produced in heavy-ion colHsions Hydrodynamics, being 

an effective macroscopic approach valid at large distances and long time scales, cannot give account of the pre- 
equilibrium stage of the reaction. However, hydrodynamics describes correctly elliptic flow at RHIC, and allows to 
jC^ ■ put some bounds on the thermalization time t < 1 fm/c. 

' An experimental evidence of the production of a thcrmalizcd quark-gluon plasma in ultrarelativistic heavy-ion 
O ■ collisions at RHIC Q is the so-called jet quenching. In this phenomenon high pT partons produced in the initial stage 
of the collision by hard scatterings loose energy and degrade, mainly by radiativeprocesses, while traveling through 
, the hot and dense medium transferring energy and momentum to the plasma (see [3| for reviews). In order to describe 
' these processes and the subsequent modification of the hadronic spectra due to the interaction of the high pT partons 
. with the medium, various models have been proposed where perturbative QCD is supplemented with medium-induced 
' parton energy loss [3l or where the AdS/CFT correspondence [H is employed. 
^ , It is interesting to analyze whether some aspects of the process of jet quenching can be described employing a 
hydrodynamical picture [6|. Hydrodynamics describes the behavior of a system in local equilibrium in terms of 
conservation laws of macroscopic quantities. On the other hand, the jet quenching is related to the modification 



^ . of high-pr parton fragmentation processes in the medium. It is then clear that hydrodynamics cannot describe the 
■ " ' microscopic processes related to jet quenching, but it can give information on the macroscopic and collective behavior 
of the system composed by the plasma and the hard jets. 

In Refs. 0, [ll the effect of a jet crossing the plasma at high speed has been analyzed. The authors have employed 
a hydrodynamical approach to describe the process of energy and momentum deposition from the fast jet to the 
surrounding plasma. In a hydrodynamical picture, a high pT jet crossing the medium at a velocity higher than the 
speed of sound forms shock waves with a Mach cone structure. Such shock waves should be detectable in the low 
Pt parton distributions at angles tt ±1.2 with respect to the direction of the trigger particle. A preliminary analysis 
of the azimuthal dihadron correlation performed by the STAR Collaboration and the PHENIX Collaboration 
(l0| seems to suggest the formation of such a conical flow. However, it is yet unclear whether those structures are 
compatible with Mach cone formation and alternative explanations have been proposed [ll[ . The study performed in 
Ref. concludes that for realistic phenomenological values of the hydrodynamical variables the Mach cone effects 
arc too weak to explain the PHENIX results. It is our aim here to explore whether other hydrodynamical effects, not 
considered in may enhance the signal. 

The study of the interaction of a relativistic stream of particles with an electromagnetic plasma is a topic of interest 
in different fields of physics, ranging from inertial confinement fusion, astrophysics and cosmology. When the particles 
of the stream are charged, plasma instabilities develop, leading to an initial stage of fast growth of the electromagnetic 
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fields. One then talks about filamcntation, two-stream or Weibel instabilities, according to which is the fastest growing 
collective mode, although the notation is non-universal and sometimes confusing. Weibel instabilities usually refer to 
the case when the plasma is out of equilibrium, and with an anisotropic distribution in the velocities of its constituents, 
but the name is also used in the context of a jet moving in an equilibrated plasma if the instability appears in the 
transverse modes. These sort of iet induced instabilities have been studied using a variety of methods, from kinetic 
theory to hydrodynamics [l^ . Il3l | . Experimental evidence of the relativistic filamcntation instability has also been 
reported in Ref. [14| . 

The study of chromo- Weibel instabilities is now a very active field of research p^-[30j (see also the recent reviews 
[sH and dl] for a more complete list of references). This is so because it was suggested [l^ US that the presence of 
plasma instabilities could be a natural explanation for the fast equilibration that measurements of the elliptic flow 
at RHIC seem to imply. In the early stage of a heavy-ion collision, the non-equilibrium anisotropic distributions of 
the partons should be responsible for the fast growth of the chromomagnetic plasma modes, which in turn would 
isotropize the system and speed up the thermalization process. Whether the experimental conditions met at RHIC 
are favoring this thermalization scenario or not is a question that requires hard numerical simulations. 

In this paper we study how a relativistic jet crossing an equilibrated quark-giuon plasma induces instabilities using 
a chromohydrodynamical approach [ssj . We will assume that both the plasma and the jet can be described using 
hydrodynamics. Therefore, regarding the plasma, we consider conditions that can be realized after t ^ 1 fm/c of the 
moment the heavy ion collision has taken place and the plasma has thermalized. Regarding the jet, there are some 
experimental evidences that the particles that constitute the jet that crosses the plasma equilibrate. Indeed, consider 
a fast parton that traverses the medium in the the direction opposite to the direction of a high px trigger particle. 
The mean value of the momentum associated of the soft hadrons emitted by the parton that crosses the medium and 
detected in the hemisphere opposite to the direction of the trigger particle reaches a common value Q suggesting that 
the energy lost by the fast parton thcrmalizcs . 

In studying the evolution of the system composed by the plasma and the jet we will employ ideal hydrodynamic 
equations that do not take into account the effect of collisions. Therefore our results will be valid on time scales 
shorter than the mean free path time for collisions. We postpone the study of more realistic fluid equations to future 
work. 

The validity of the chromohydrodynamical approach will be our starting assumption, as this allows us to simplify 
the equations governing the evolution of the system. To the best of our knowledge, only Ref. fs^l considers the 
possibility of the appearance of fllamcntation instabilities produced by hard jets in heavy-ion collisions. However, the 
approach considered there and ours are different. The analysis of Ref. [s^l is performed within kinetic theory, and 
thus relies on the quasi-particle picture and a weak coupling scenario. Instead, we use fluid equations. We note that 
colored plasma waves are believed to be very quickly damped (and this is the reason to exclude them in the study 
of Ref. (6[). That conclusion might have to be reviewed, as in the presence of instabilities rather than damping, one 
finds exponential growth of the fields. Our approach is also different to the color wake field scenario of Refs. [H, [11]. 
In those references, the authors studied the field response to a single moving colored particle, which produces a color 
charge density wake. Instead, we consider a color neutral jet, treated hydrodynamically, that moves trough the plasma 
and produces color fluctuations and instabilities. 

This paper is structured as follows. In Section |IT] we review the chromohydrodynamical equations describing 
the fluctuations of a plasma around the stationary colorless state in the linear approximation [33| . In Section IIIII 
the equation describing sound colorless fluctuations and plasma colored fluctuations are derived. In Section IIVI we 
consider the case where a relativistic jet crosses the plasma. Assuming that both the plasma and the jet can be treated 
using a hydrodynamical approach, we derive the dispersion law for the collective modes. We find that there is one 
unstable mode if the velocity of the jet is larger than the speed of sound and if the momentum of the collective mode 
is in modulus smaller than a threshold value. Regarding the orientation of the momentum of the collective mode, we 
study separately the case where it is coUinear with the velocity of the jet, the case where it is orthogonal and then 
for a generic angle between the two. Quite interestingly we find that the unstable modes with momentum parallel to 
the velocity of the jet is the dominant one for velocity of the jet v < 0.8. For larger values of the jet velocity only the 
modes with angles larger than ~ tt/S are significant and the dominant unstable modes correspond to angles ^ 7r/4. 
We draw our conclusions in Section [V] We work using natural units h = c = ks = 1 and use metric convention 
(1,-1,-1,-1). 



II. CHROMOHYDRODYNAMIC EQUATIONS FOR THE QUARK-GLUON PLASMA 

Hydrodynamical equations are the expressions of the conservation laws of a system when it is in local equilibrium. 
In the quark-gluon plasma there are conservation laws which concern the baryon current, the color current and the 
energy-momentum tensor. In Ref. [sj the local equilibrium state for the quark-gluon plasma has been determined. 
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It is in general described by one singlet four velocity, a baryon density, singlet energy and pressure, and in principle, 
a non-vanishing color density. However, dynamical processes associated to the existence of Ohmic currents tend to 
whiten the plasma quickly, on time scales much shorter than momentum equilibration processes (38| . For this reason 
one can expect that only colorless (singlet) fluctuations are relevant at large time and space scales. However, there 
are situations, as the one considered in the present paper, when color fluctuations grow on short time scales instead 
of being damped. Therefore in order to describe the short time evolution of the plasma, one needs to include color 
hydrodynamical fluctuations in the equations. In Ref. [33j such a chromohydrodynamical approach for the short time 
evolution of the system has been formulated, and here we will review the basic set of equations. 

For simplicity, as in Ref. [ssj , we will only consider the contribution of quarks in the fundamental representation. 
The inclusion of antiquarks and gluons is straightforward. The fluid approach is based on the covariant continuity 
equation for the fluid four-flow 

D^n^^ = (1) 

and on the equation that couples the energy-momentum tensor T^^ to the gauge fields 

i?pT^''-|{F/,n^} = 0, (2) 

where the various quantities are 3x3 hermitian matrices in color space, and we have suppressed color indices. The 
covariant derivative is defined as 

Df, = d^- ig[Afj,{x), ... ] , 

with Af^ ~ A^{x)t°' for a ~ 1,...,8 and t'^ — A°/2, where A" are the SU{3) Gcll-Mann matrices and therefore 
Tr(T°T^) = |<5"''. The strength tensor appearing in Eq.Q is given by F^^,^ = df^A^ — d^A^ — ig[A^, A„]. 

We further assume that the four-flow and the energy-momentum tensor have the expression valid for an ideal fluid, 
i.e. 

n^'{x) = n{x)u^'{x) , (3) 

and 

T^'i^) = li^i^) +Pix)){u'^ix),u^x)}~pix)g^'' , (4) 

where the hydrodynamic velocity w^, the particle density n, the energy density e and the pressure p are 3x3 matrices 
in color space. The quantities defined above have in general both colorless and colored components, as an example 
the particle density can be written as 

nap = riolafj + -riaT^p , (5) 

where a,/3 = 1,2,3 are color indices and / is the identity matrix. In the following equations we will omit the color 
indices not to overcharge the notation. 

The color current due to the flow of the fluid can be expressed in terms of the hydrodynamic velocity and the 
particle density as 

f{x) = -l^(nu^-h:r[nu^]), (6) 

and it acts as a source term for the gauge fields in the Yang Mills equation 

D^F>^^{x)^f{x). (7) 

Thus, we will assume that all the gauge fields that appear in the fluid equations are only due to the presence of a 
colored current in the medium. 

It is worth remarking that Eqs. ([T|) and ^ were derived in Ref. [11] from the coUisionlcss transport equation obeyed 
by the particle distribution function. Thus, the conservation laws expressed by the Eqs. ([T]) and ^ are strictly valid 
on time scales shorter than the mean free path time. As an example, for quarks rP represents the quark particle 
density, which fulfills a conservation law that says that the change of the number of particles within a volume clement 
is equal to the fiux of particles across the surface of the volume element. On the other hand, for times larger than the 
mean free path time this conservation law is violated by collision processes that change the particle number inside 
the volume element. Then, only baryon number is conserved, which requires the knowledge of both the quark and 
antiquark particle densities. 

Summarizing, for time scales shorter than the mean free path time there are more conservation laws that at long 
time scales, as Eqs. ([Ij and ^ indicate. Their validity for the short time phenomena that we will be studied here is 
then guaranteed. 
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A. Linearization of the chromohydrodynamic equations 



We shall consider the fluctuations of density, energy density, pressure and plasma velocity, around their stationary 
and colorless state described by n, e, p and u'^ respectively. In order to study such fluctuations we linearize the 
chromohydrodynamic equations ([Tl[2), assuming the fluctuations to be small. Notice that in the stationary state the 
color current defined in Eq.® vanishes, indeed 

jf" = -g (^n u'" - ^Tr [n u^'] ^ = > (8) 

and we will assume that in the stationary state no field F^'' is present in the system. We define the space dependent 
fiuctuations of the various quantities around their colorless values as 

n{x) = n + Sn{x) , e{x) = e + Se{x) , (9) 



p{x) ^ p + Sp{x) , u^'{x) ^uf" +Su''{x) . (10) 
All the fluctuations can contain both colorless and colored components, therefore they can be decomposed as 

Suap = Snolafj + ^SUaT^p , (11) 

where Sno is the colorless fluctuation, and 6na the corresponding colored components. 

The state described by n, e, p and is assumed to be stationary, colorless, and homogeneous on the scale of 
variation of the fluctuations and therefore we have that 

Df'n = , D^e = , D^p = , Df^u" = . (12) 

Moreover, since we will consider only small deviations from the stationary state, we will assume that the following 
conditions are obeyed 

n > (5n , e > 5e , p:$> Sp , u'' > Su^" . (13) 

Actually, Sn, Se, Sp and Su^^ should be diagonalizcd to be comparable to the n, e, p and u'^. 

We now aim to determining the set of equations that the fluctuations Sn, Se, Sp and Su'^ obey. Employing the 
equations ^ and PU)) we can derive the expression of the fluctuation of and of T^'' around their colorless and 
homogeneous values: 

n'^ =nu'^ + n5u'^ +6nu'' , (14) 



T^"' ^ {€+p)u^'u'' -pg^'•' + {S€ + Sp)u^'u'' + {€+p){u^'Sll'' + Su^'u'')-Spg^'•' . (15) 

Upon using these expressions in Eqs. mH]), we obtain that Su^^, Sp, Se and Sn must obey the following equations: 

nDf.Su^' + {D^Sn)u^' = , 
u^'Df,Se + {e+p)D^M'" = 0, (16) 
{e+p)Uf,D>'Su'' - {D" ~ u''ui,D'')Sp- gnUi.F'"' = 0. 

These equations do not form a closed set and one more relation has to be provided. In hydrodynamical treatments, 
one usually imposes an Equation of State (EoS) of the form 

p = p{e,n), (17) 

and in this way one obtains a closed set of equations. For the applications wc have in mind, we will consider that 
the conformal limit is reached, and further the effects of the particle density can be ignored. Thus wc will use an 
Equation of State given by 

p{x)=cleix), (18) 

where Cg is the speed of sound, and in the conformal limit ~ l/\/3- We will however leave as a parameter, and 
use its value only in the numerical analysis of the equations. 
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We will study separately colorless and colored fluctuations. In order to obtain the equation governing the colorless 
fluctuation we take the trace of the Eqs. obtaining 

ndfj,5u'^ + {d^Sno) u'' = , 
u'^d^Seo + (e + p)d^Su^ = , (19) 
{e + p)u^d^Su'^ - id" - ru^d^Spo = . 

Whereas multiplying the Eqs. (jl6p by r" and taking the trace we gather the equations for colored hydrodynamical 
fluctuations: 

u''{D^Se)a + {l + p){D^Su^)a^O. (20) 
{e + p)u^iD^5una - {{D- - u-u^D^^)5p)a - gnu^F^ = , 

where we have used the notation {Df,X)a = d^X'' + gf^'^A'^^X''. 

III. SOUND AND PLASMA WAVES IN THE QUARK-GLUON PLASMA 

In order to analyze the linearized chromohydrodynamical equations derived in the previous Section we will treat 
separately the colorless and colored fluctuations showing that they describe the propagation of sound and plasma 
waves respectively. 

A. Sound waves 

The colorless fluctuations of plasma velocity, energy density, pressure and density obey the Eqs. (fT9|) that in 
momentum space read 

u^k^Sno + nk^SuQ ~ , 
u'^fc^(5eo + (e + p)A.v(5< = 0, (21) 
i {e + p) Up^k'^SuQ + i{u'^u'^k^ ~ k^)6pQ = 0, 

where fc'^ = (cj, k). Since pressure and energy density are related by the EoS (|18p . their variations are not independent 
and upon differentiating Eq. (|18p we obtain that 

Spo^ (^)seo = cl5eo . (22) 



Substituting this expression in Eqs. (|2ip we find the following equation for the fluctuations of the pressure 

1 



u'^k^ 1 \ c? 



^ - 1 {u^k^y + k'}6pa^0, (23) 



while the fluctuations of the other quantities can be expressed as 



5no = -nM^. (25) 
u ■ k 



In the plasma rest frame, u'^ = (1, 0, 0, 0), Eq. simplifies to 



1 k' 



tj2 



Spo = 0, (26) 



which gives the standard expression of the sound waves in a plasma, reflecting the fact that all colorless hydrodynamical 
fluctuations propagate at the speed of sound. 
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B. Plasma waves 



Regarding the colored components of the fluctuations, Eqs. ((20)) . we hnearize in the A flelds as wch. The reason 
is that we will be interested in the computation of the polarization tensor and expanding at the first order in A is 
sufficient. One then finds 



u^'kf,Sea + ie+p)kf,Sui^ = 0, 
t{e+p)u,,k^Su''^+i{u''u''k^-k'')Spa-qnu,,F^'' ^ 0, 



(27) 
(28) 
(29) 



and the fluctuation of the pressure and the fluctuations of the energy density can be related using the EoS (fT8|) , that 
leads to 6pa = c^Sea p^,[39^. 

Equations (P7)) and (^5)) relate the fluctuations of the density and energy density with the fluctuations of the velocity 



Sua 



u ■ k 



and combining these expression with Eq. ([22]) we find that Eq. ((29|) can be rewritten as 



(30) 
(31) 



:{k''k'' -u-^k^iu-k)) 



5u,, 



{e + p){u-k) ^ 



{u-kf 

The inverse of the operator on the left hand side of this equation can be determined observing that 



{u ■ ky 



where 



(fc'^fc^ - u'^k'^iu ■ k)) - T{k){k^k^ - u„k^{u ■ k)) 



T{k)^~- 



F + (^-l)(M-fc)2 ■ 

Then one can solve Eq. ([5^ and the fluctuation of the hydrodynamic velocity takes the form 



<5u<T,a = ig 



{e + p){u ■ k) 



g<yv + T{k) {kaku — Uakv{u ■ k)) 



(32) 



(33) 



(34) 



(35) 



The fluctuation of the current induced by the fluctuation of the density and of the hydrodynamic velocity can be 
derived from Eq.fG]) and is given by 



Such fluctuations are related in linear response theory to fluctuations of the gauge flelds via 

Sj^ik) = -n;:afc)A,,,(fc) , 



(36) 



(37) 



where H^^ is the polarization tensor. Considering that the linearized strength tensor equals Fj^'^{k) = —ik'^A'^{k) + 
ik^Ai;^{k), and upon substituting the values of the fluctuations of density and energy density in Eq.([36]) we obtain 
that 



1 



2 (e + p) {u-ky 



{u ■ k){k>'v!' + k^u^") - {u ■ k^g'"' - k'^u^u" 
1 



((iZ • k)k'^{k''u'' + k^'u") - {u ■ k^k^k" - k^f'u''^ (.38) 



fc2 + (-|-l)(i2.fc)2 

One can easily check that the polarization tensor ([55)1 is symmetric and transverse (fc^n^'^(fc) = 0). 
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In the plasma rest frame, = (1,0,0,0), the polarization tensor simplifies considerably and its components are 
given by 



fc2 



k2 



r,2 ij U 
f\j fx 



C2k2 ' 



(39) 



C2k2 



where 



2 (e+p) 



(40) 



is the plasma frequency. 

At the linear order, the Fourier transformed chromodynamic field j4^(fc) satisfies the equation of motion 

fc2g*'''-A;^r-n^''(fc)jA^(fc) = 0, (41) 

where we have dropped the color indices a, b. 

In order to determine the dispersion law for the gauge fields we define the dielectric tensor e*-' (fc) according to 



ey(fc) = 5^J + ^ff^(A:), 



that in the plasma rest frame is given by 

4'(^'k) 



1 



up- up- — c^k^ 

The dispersion equation can now be determined by computing the determinant of the matrix 



M'^ = - k'k^ - uj'e'^ (fc) , 



and solving the equation 



det 



0. 



(42) 



(43) 



(44) 



(45) 



where we assumed the Coulomb gauge, k • A{k) ~ 0, and ~ 0, meaning that E = iujA. 

Taking the wave vector in a given direction, say k = (fc,0, 0), one finds the following solutions to the dispersion 
relation 



,2 
P 



2 7,2 

CgK , 



(46) 
(47) 



which correspond to longitudinal and transverse modes, respectively. Notice also that both modes are not damped. 
The reason is that we have neglected the dissipative terms in the hydrodynamical equations, assuming that the colored 
plasma on the time scale of interest to the present analysis can be approximated as an ideal fluid. In order to consider 
the damping of color fluctuations in the hydrodynamical limit one should include transport coefficients, such as color 
conductivity or color diffusion. 

Let us compare the dispersion laws of the longitudinal and transverse modes obtained above with the corresponding 
results in kinetic theory. Using the hard thermal loop (HTL) effective theory one can determine the polarization tensor 
for the gauge fields. The dispersion laws in the limit fc <C for the longitudinal and transverse modes are given by 

m 



(48) 
(49) 
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respectively. The chromohydrodynamical approach reproduces the short time physics of kinetic theory in the domain 
k <^ Lu. There is a numerical discrepancy in the terms when comparing the two set of dispersion laws. The origin 
of such a discrepancy is that in order to go from transport theory to a fluid approach, one has to impose a relation 
between the energy and the pressure of the system that allows to discard all the higher momenta moments that are 
contained in the particle distribution function. Such an approximation, typical of an hydrodynamical approach, does 
not affect in any profound manner the physics described in the domain of applicability of the formalism. 



IV. A RELATIVISTIC JET TRAVERSING THE PLASMA 



We shall now consider a jet of particles propagating across the plasma at a relativistic velocity v — |v|. We assume 
that both the jet and the plasma are initially in equilibrium but jet and plasma are not in equilibrium with each other. 
The system formed by the plasma and the jet is one specific example of the so-called two stream systems, which have 
been widely studied in magnetohydrodynamics. and extended to the chromohydrodynamical case in Ref. [33| . 

We are not concerned with the colorless hydrodynamical fluctuations, which have been studied elsewhere 0, H]. 
We will instead examine the evolution of the colored fluctuations at short time scales, where the equations describing 
such fluctuations can be linearized. 

We consider the frame where the plasma is at rest. The velocity, as well as the pressure, energy density and particle 
density of the jet and of the plasma will be different and will be labeled in a different way. Since the plasma is at rest, 
uj^ = (1,0,0,0); whereas the velocity of the jet will have the general expression uj^j = 7(1, v), with 7 = 1/ \/l — v^. 
The polarization tensor associated to the jet can be computed in an analogous way to that of the plasma and will 
have the expression reported in Eq. (|38|1 with -S^ — > -Uj^j. In the limit where 7 3> 1, it is possible to neglect the second 
term in the square bracket on the right hand side of Eq. ((38|) . Such an approximation corresponds to solving the fluid 
equations neglecting the effect of the pressure gradients. One can also see that this approximation corresponds to 
considering that the distribution function of the constituents of the jet is of the form (see [33| for a discussion of this 
appr oximat ion ) 

f(p)^nuU'-'^{p-'-±Pn) , (50) 

which describes fast moving particles, with a non-thermal distribution. 

Employing Eq. we obtain the following expression of the dielectric tensor for the jet: 



. (1 - ^)s" _ ^ii(^™! _ (f-f'-'f 

V Lu^ J Lu^ \ uj-k-v (w - k • v)2 



(51) 



where 



2 „ 5' 



2 (Cjet+Pjet) 

is the plasma frequency squared of the jet. 

In order to obtain the dielectric tensor of the system composed by the plasma and the jet we note that the total 
polarization tensor of the system is given by the sum of the two polarization tensors in linear response theory 

nr-ni^'^ + nf,';, (53) 



where the expression of the components of TM^" is given in Eq. ((39)) . The dielectric tensor of the total system turns 
out to be 

where 



Lo"^ I uP- llP- — Cgk^ \ cij — k • V (u; — k • v) 



'A = + ^fet ■ (55) 

The dispersion laws of the collective modes of the system composed by the plasma and the jet can now be determined 
solving the equation 



det Ml^ = det 



^^8'^ -k'¥ -u?e^{k) =0. (56) 
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The solutions of this equation depend on |k|, |v|, cos9 = k • v, Wp, Ujet and Cg. We find that Eq. (|56p admits eight 
solutions, however two of them are trivially given by w = and six of these solutions are given by the roots of a 
polynomial of the sixth order. In order to simplify the expression of the dispersion law and to study the dependence 
of the unstable collective modes on the values of the various parameters, we define the dimensionless quantities 

y--, 5=^%, (57) 

Wt UJt Lot 

meaning that we will measure momenta and energies in units of Ut- Since the plasma frequency of the jet is unknown 
we will treat 6 as a parameter and we will analyze values of 6 ^ 1 corresponding to LOj^t Wp. 

We anticipate that we find that five of the dispersion laws are always stable, whereas one is unstable in a certain 
range of values of the parameters. In particular we find that for values of v smaller than the speed of sound of the 
plasma, Cs, all the modes are stable independent of the values of the remaining parameters. Let us also note that 
when k and v are not parallel, it is not possible to decompose the dielectric function (|54p with only transverse and 
longitudinal projectors of k. We will then not refer to transverse or longitudinal modes in that case. 

In the following Subsections, we will analyze the dispersion relations of the collective modes for different orientations 
of k and v. 



A. k parallel to v 

Here we consider the case corresponding to k || v, and choose k = (0, 0, k) and v = (0, 0, v). This case can be easily 
analyzed because for this orientation of k and v the dielectric tensor ([M)) and the matrix M^-' defined in Eq. ([5^ 
are diagonal and the corresponding dispersion equation factorizes. Moreover, since k is parallel to v it is possible 
to define transverse and longitudinal modes with respect to the orientation of these vectors. We find two stable 
transverse modes with dispersion law 

Lo^ = ujI + fc2 . (58) 

Notice that this dispersion law is analogous to the dispersion law that we obtained in Eq. ([T7|) that describes the 
propagation of the transverse mode in a a plasma without a jet. The only difference is that the plasma frequency 
Lo-p has been replaced with Wf Therefore the only effect on the transverse modes due to the presence of the jet is to 
change the plasma frequency of the mode. 

Regarding the longitudinal modes, they are given by the solution of the following equation: 

2 2 2 c^^^ 2 '2'Vkw — v'^uj'^ — v'^k'^ 

- - - ^T^H^ = ° • 

In terms of the dimensionless variables defined in Eq. ([57)) . the equation for the longitudinal modes (|59p can be written 



as 



- 2vyx^ + (-1 + bv^ + 2/2(„2 _ c^)) + 2vyx{l -b + cly^) - {bv^icl - 1) - c^fe + v^ic^y^ + 1)) = 0. (60) 



This equation has two trivial solutions x ~ 0, and four solutions corresponding to the roots of a quartic equation. As 
already mentioned we find that at most one mode is unstable. In the present case one of the longitudinal modes is 
unstable for Cs < w < 1. Since the cases where v = Cs and w — > 1 will play a special role for these modes, let us first 
consider the solutions of Eq. ([5^ for these two values of the velocity. 
In the limit v ^ 1, Eq. ((59)l simplifies to 



with solution, cj^ = and 



„2i,2 2jL2 
'^t P^2_c2fc2 ^'^jct ^ p P^2_c2fc2 ! 



uj' =uj'p + cik' . (62) 

This solution is the same that we have obtained for a plasma without a jet in Eq. (j46p . Therefore in the limit v ^ 1, 
the longitudinal mode is stable and we see that a jet of particles moving at the speed of light can only affect the 
transverse modes , but they do not affect at all the longitudinal modes. This is related to the dimensional contraction 
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which occurs in the eikonal Hmit [4l| , when it can be shown that the strength field tensor caused by a massless particle 
that moves at the speed of light is on the plane transverse to its motion. 

In the case v = Cg it is possible to find one simple analytical solution (apart from the trivial solution uj^ ~ and 
Eq. ([55]) ) which is given by 

uj = Cgk , (63) 

whereas three more solutions correspond to the roots of 

- CsX^y + Csvil + b{-2 + c^) + a;(-l + cl{b - y^)) + cly^) = . (64) 

The discriminant of this equation is negative for 6 < 1 meaning that the three solutions of this equation arc always 
real and the corresponding modes stable. 

We now turn to a numerical study of the non-trivial solutions of Eq. (|60p as a function of the parameters b and v. 
We find a solution with a positive imaginary part for any non- vanishing value of h and for v larger than Cs- Such a 
solution corresponds to the unstable longitudinal collective mode of the system. In Fig. [1] we have reported the plot 
of the imaginary part of the frequency, F, of the this longitudinal mode as a function of k/ujt^ for h = 0.02 and for 
four different values of the velocity of the jet: v = 0.6 full (red online) line; v ~ 0.7, dashed (magenta online) line; 
V = 0.8, full (green online) line; v = 0.9, dot-dashed (blue online) line. For any value of the velocity smaller or equal 
the speed of sound we numerically find F = 0. For values of v larger than the speed of sound and less than 1 this 
collective mode becomes unstable in a range of values of the momentum k < fcmax- For small values of fc, F increases 
with increasing momentum and reaches a peak. Then it decreases and at fcmax becomes zero. As can be seen in Fig. [1] 
for V > Cs, the value of fcmax decreases with increasing velocity. 

The value of F at the peak, that we will indicate with Fmax, depends on v and b. In Fig. [5] we present the plot of 
Emax as a function of b for four values of the jet velocity: v = 0.6 full (red online) line; v = 0.7, dashed (magenta 
online) line; v = 0.8, full (green online) line; v = 0.9, dot-dashed (blue online) line. Note that independent of the 
value of V, for 6 = 0, we obtain that F = 0. In this case the plasma frequency of the jet is zero and corresponds to 
the case where the system consists of the plasma only. Indeed, for & = 0, the dielectric tensor ((54|) does not have any 
contribution from the jet and therefore all the collective modes are stable, as shown in the previous Section. With 
increasing values of 6, ie. of the plasma frequency of the jet, F^ax increases, meaning that to larger values of the 
density of the jet (for fixed values of g and ejot +Pjct) correspond larger instabilities. Moreover, for values of & < 0.01 
the value of Fmax increases quickly with increasing b whereas for larger values of 6, Fmax becomes less sensitive to the 
actual value of h. In particular for v = 0.6, F^ax saturates at ^ 0.09 As a function of the velocity we find that 
for a given value of b the value of Fmax increases with increasing velocity for Cs < v < 0.7. For v ~ 0.7 — 0.8 the 
value of Fmax reaches its maximum value and then for larger values of the velocity decreases and eventually in the 
ultrarelativistic case v — ^ 1 we find that the frequency of the longitudinal mode becomes real for all the values of k. 



— v = 0.6 




FIG. 1: (Color online) Imaginary part of the dispersion law of the unstable longitudinal mode for the system composed by 
a plasma and a jet as a function of k/ujt- Here the momentum of the unstable mode is parallel to the velocity of the jet v, 
h = ojj^ct/'^t = 0.02 and the four different lines correspond to different values of the velocity of the jet |v|. 
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FIG. 2: (Color online) Variations of the largest value of the imaginary part of the dispersion law of the unstable longitudinal 
mode for the system composed by a plasma and a jet in the case k || v as a function of b for four different values of the velocity 
of the jet, |v|. 



B. k orthogonal to v 

We now consider the case corresponding to k orthogonal to v, choosing k = (fc, 0, 0) and v = (0, 0, v). In this case 
the matrix M^'' defined in Eq. (|56p is block diagonal. One block is one-dimensional and corresponds to the component 
Mf^. The other block is two-dimensional and corresponds to the components of AI^-' with indices i,j = x, z. Therefore 
the dispersion equation (|56p factorizes into two equations. Regarding the one-dimensional block it describes a mode 
orthogonal to both k and v. This mode is stable and has dispersion law 

Lu^ = uj^ + . (65) 

This dispersion law is analogous to the dispersion law for the transverse mode for a plasma without a jet, given in 
Ea. (|47p . with the plasma frequency LUp replaced with ut- The remaining collective modes are solutions of the equation 

Using the variables defined in Eq. (j57p . one can rewrite this equation as 

{x^-1- cly^) (x^ - 1 - j/2 + ^,2^) + x^y%{cl - v^) -f hy^ [v^ -cl + cly^{v^ ~ 1) + clv^b) - b^v^y^ = , (67) 

where we have already factored out the trivial solution x'^ = 0. 

Also in this case we find that for v < Cg no mode is unstable and for v > Cs only one mode is unstable in a certain 
range of values of the momentum. We have numerically studied the solutions of Eq (j67p and we have reported the 
results in Fig. [3] and Fig. 21 In Fig. [3] it is shown the plot of the imaginary part of w for b = 0.02, as a function of k for 
five different values of the jet velocity: v ~ 0.6, full (red online) line; v = 0.7, dashed (magenta online) line; v ~ 0.8, 
full (green online) line; v = 0.9, dot-dashed (blue online) line; f = 1 dot-dot-dashed (black) line. 

The largest value of k where the orthogonal mode is unstable is a monotonic increasing function of v and it diverges 
for w — > 1. Also the value of the imaginary part of F at the peak, Fmax, increase with increasing v, but for v ^ 1 
reaches a maximum finite value. The behavior of the largest value of the imaginary part of the frequency, Fmax, as a 
function of 6, for various values of the velocity, is reported in Fig. 3) With increasing values of 6, i.e. with increasing 
value of the plasma frequency of the jet, the value of r„iax increases. Moreover, for a given value of b, the larger the 
velocity of the jet. the larger is the value of Fmax- 



C. Arbitrary angles 



In the previous two Subsections we have analyzed the behavior of the unstable mode of the plasma for two values 
of the angle 9 between k and v. In particular in Section IIV Al we have considered = 0, corresponding to modes 
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FIG. 3: (Color online) Imaginary part of the dispersion law of the unstable mode for the system composed by a plasma and a 
jet in the case where the momentum of the collective mode, k, is orthogonal to the velocity of the jet, v, as a function of the 



momentum k/ujt and for b - 



/ujt = 0.02. The five different lines correspond to different values of the velocity |v| 
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FIG. 4: (Color online) Largest value of the imaginary part of the dispersion law of the unstable mode for the system composed 
by a plasma and a jet in the case where the momentum of the collective mode, k, is orthogonal to the velocity of the jet, v, as 
a function of 6 = Lu^^^jiJl. The five different lines correspond to different values of the velocity |v|. 



with k 11 V and in Section llV Bl we have taken Q = 7r/2, corresponding to modes with k _L v. A common feature of 
the two cases is that the value of Fmax (corresponding to the largest value of the imaginary part of the frequency as 
a function of |k|), increases with increasing values of the plasma frequency of the jet. Regarding the dependence of 
Tmax on the value of the velocity of the jet v one can see, comparing the plots in Fig. [4] with those in Fig. [2l that 
for ultrarelativistic jet the modes with k _L v dominate the modes with k || v. On the other hand, for velocity larger 
than the speed of sound, but v < 0.9 the longitudinal modes are dominant. 

In this Subsection we analyze the general case, of an arbitrary angle d and we will determine at which angle 
corresponds the most unstable mode as a function of v and h. For definitiveness we choose v = (0, 0, v) and k = 
(0, fcsin6', fccosfl). 

The analysis of the dispersion laws for the modes with k || v and k _L v was simplified by the fact that in both 
cases the matrix defined in Eq. ijSS]) becomes block-diagonal and the corresponding equation for the dispersion 
laws factorizes. However, for arbitrary angle d the matrix is not block-diagonal. Then we will rely on a numerical 
solution of Eq. (j56p . Employing the definitions in Ea. (|57p we obtain (apart from the two trivial solutions corresponding 
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to = 0), that the six dispersion laws can be obtained from the roots of the equation 

= 2x^ + x^{-A + 26w^ + (-2 - + v^)y^) +x^{2~ 2bv^ - (-2 - 2(1 + b)cl + {2 + b + 2bcl)v^)y'^ 
-{v' + cl{^2 + 1,2))/) + by\~{v\b + y2)) + ^^(-2(1 + y^) + v\l + b + 2y^))) 
+vy{v{l + cly'^){y + y^) - 4x(-l + ~ y'^){~l + b + x^ - cly'^) cos{9) 

+vy{-{b''{^l + cD) + b{-2 + cl + x^- y2) + (_l + ^2 _ y2^(^_^ ^ ^2 „ ^2^2)) cos(20)) . (68) 

Of the six solutions of this equation only one corresponds to an unstable mode in a certain range of parameters. In 
particular, for any value of 6 wc find an unstable mode only if v > Cg and if the momentum is smaller than a certain 
value fcmax, where fcmax depends on 0, v and b. The largest value of the imaginary part Fmax depends on 9, v and b as 
well and in Fig. [5] we show three plots of Fmax as a function of b for various values of 6 and v. In each plot the full (red 
online) line corresponds to 9 — 0, the dashed (magenta online) line corresponds to 9 ~ tt/8, the dot-dashed (green 
online) line corresponds to 6* = 7r/4, the dotted (blue online) line to 6* = Btt/S and the dot-dashed-dashed (black) line 
to 9 = Tx 12. The left panel corresponds to u = 0.8, the central panel to w = 0.9 and the right panel corresponds to 
w = 1. 

For V = 0.8 (left panel), the most unstable modes are those corresponding to small angles 0, i.e. are those modes 
with k almost coUinear with the velocity of the jet. In this regime these modes dominate the dynamics. With 
increasing velocity modes with larger values of the angle become relevant. Indeed for a value of the velocity v = 0.9, 
corresponding to the central panel in Fig. O the imaginary part of the modes with < < 7r/4 is the largest and 
the corresponding modes dominant. In the ultrarelativistic case t; — !■ 1, right panel in Fig. [5l the dominant modes 
are those with 9 ~ 7r/4. Modes with k almost coUinear with v, or more precisely all the modes at angles 9 < 7r/8, 
are suppressed. Notice that the modes with very large values of ~ 7r/2 are not the dominant one. This means that 
modes with k orthogonal to v are not the most important: more important are "oblique" modes with ~ 7r/4. 

In agreement with the analytical results of the previous Section, in the limit w — + 1, that is shown on the right panel 
of Fig. [51 the mode with 9 = Q has a vanishing imaginary part and is stable. 




FIG. 5: (Color online) Largest value of the imaginary part of the dispersion law for the unstable mode as a function of 6 for 
three different values of the velocity of the jet |v| and five different angles between k and v. The left panel corresponds to 
V = 0.8 and the dominant unstable modes correspond to small angles 8, i.e. to modes that are almost coUinear with the jet. 
The central panel corresponds to v — 0.9 and the right panel to v —> 1. In both cases the modes corresponding to S ~ 7r/4 
are the dominant one. In the limit v —> 1 modes with momentum coUinear with the velocity of the jet are suppressed and the 
mode with 6' = vanishes. 



V. CONCLUSIONS 



It is well-known that hydrodynamics predicts that any object moving in a fluid at a velocity higher than the speed 
of sound creates shock waves and a Mach cone structure. The conclusion of our study is that a neutral stream of 
colored particles moving at a velocity higher than the speed of sound in a quark-giuon plasma in the conformal limit 
also causes gauge field instabilities. It is curious to note that it is the speed of sound the threshold value for the 
appearance of such instabilities. Qualitatively one can understand why this happens considering that if an external 
perturbation to the plasma moves slower than the speed of sound, the system can respond by locally rearranging 
the values of density, energy density, pressure and plasma velocity, and then hydrodynamical fluctuations of these 
quantities may counteract the slow external perturbation. This is not the case when the perturbation propagates 
faster than the larger velocity of propagation of hydrodynamical fluctuations, i.e. larger than the speed of sound. 
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Even if we have obtained our results by solving the chromohydrodynamical equation of the system in the linear 
response approximation let us stress that the same results are easily translated for electromagnetic relativistic fluids. 
To the best of our knowledge, the speed of sound bound that we found has not been discussed in the existing literature. 

We believe that the mechanism proposed in this article provides one possible effective hydrodynamical description 
of the jet quenching phenomena, valid at macroscopic scales. Indeed we have shown that the energy and momenta 
stored in the total system (composed by the plasma and the jet) is effectively converted into (growing) energy and 
momenta stored in the gauge fields, which are initially absent. The jet looses energy in its dynamical evolution. 

However, our analysis of this mechanism has to be considered as preliminary and therefore wc cannot indicate 
consequences for existing data of jet quenching at RHIC. Indeed, we have neglected several aspects of the system 
produced in heavy-ion collisions that complicate the treatment of the jet quenching phenomenon, such as the expansion 
of the plasma, or the transition to an hadronic phase, that takes place when the plasma becomes sufficiently dilute. 
One of the relevant points here is whether the system has enough time to generate the growth of the gauge fields before 
hadronization begins. As found in our numerical study, the maximal value of the growth rate, here characterized by 
Tmax, runs between Fmax ^ (0.08 — 0.15) cjt. Therefore, the plasma instabilities fully develop on time scales of the 
order t ^ (6.7 — 12.5) /iut- To get an upper bound of that time scale, we evaluate the plasma frequency in a weakly 
coupling scenario at T ~ 350 MeV, finding t ~ 1 — 2 fm/c. 

Another aspect that wc have neglected is the existence of dissipativc terms in the hydrodynamical equations that 
at a certain stage of the plasma evolution might be able to damp the hydrodynamical fluctuations. 

In any case, assuming that the generated gauge fields have sufficiently time to grow, what is their fate? How do they 
hadronize? The late stage evolution of the chromohydrodynamical instabilities may depend on a different number of 
facts. One should discuss what is the saturation mechanism of the instabilities, which could be drastically affected 
by dissipative color damping phenomena. In any case, naively speaking, we would expect that the existence of some 
gauge field modes at soft values of the momenta less than fcmax, and in a given direction in momenta space, should 
be translated at RHIC into an enhanced production of hadrons with momenta < fcmax, versus the case without gauge 
field instabilities. Moreover, if one takes into account the effect of the growth of the gauge fields in the numerical 
simulations of the conical fiow then this might lead to an enhancement of the number of soft hadrons produced by 
the jet with respect to the case where the effect of the growth of the gauge fields has not been considered 0. 
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